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Abstract 



In this paper we solve for the quantum propagator of a general time de- 
pendent system quadratic in both position and momentum, linearly coupled 



o 

<N 
O 
CO . 

, to an infinite bath of harmonic oscillators. We work in the regime where 

q-i the quantum optical master equation is valid. We map this master equation 

(3JT), to a Schroedinger equation on Super-Hilbert space and utilize Lie Algebraic 



techniques to solve for the dynamics in this space. We then map back to 



the original Hilbert space to obtain the solution of the quantum dynamics. 
The Lie Algebraic techniques used are preferable to the standard Wei-Norman 
methods in that only coupled systems of first order ordinary differential equa- 
tions and purely algebraic equations need only be solved. We look at two 
examples. 
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I. INTRODUCTION 



In this note we describe the solution of the quantum optical master equation for a general 
time dependent system at most quadratic in position and momentum. The coupling to the 
bath is taken to be a linear and the bath is comprised of harmonic oscillators. In the optical 
regime we can write an approximate quantum master equation H 

p=-~[H,p]+Ap , (1.1) 

Ap = 7( ^ 2 + ^ {Kpa f ] + [ap,a f ]} + y {[« f ,H + Wp,a]} , (1.2) 

where 7 is the strength of the coupling to the bath and n = [exp(— (3huj) — 1] _1 is the equi- 
librium number of photons in the bath. The approximations necessary to obtain equation 
( |1.1|) are: rotating wave approximation (RWA), small system anharmonicity, Markovian 
evolution, weak coupling and initially no correlations between the system and bath. It is 
generally accepted that equation fllTTl ) describes well the physics near optical frequencies and 
forms the basis of many models studied in quantum optics [|l]|| (and references therein). 
We take as our system the Hamiltonian 

H/h = u{t) a^a + /i(t) a + a f + f 2 (t) a 2 + f 2 (t) a t2 , (1.3) 

where u, fi, f 2 , fx, J 2 are arbitrary functions of time. We will assume to G R but will not 
necessarily assume that /* = f\ and f 2 = f 2 - This is to accommodate the phenomenological 
modeling of dissipation through a non-unitary Hamiltonian 0. The general form of the 
system ( |1.3|) covers many specific models including the simple harmonic oscillator, the para- 
metric oscillator, the forced harmonic oscillator and the degenerate parametric oscillator. 
The isolated system also exhibits many families of coherent states (a la Perelomov |3J]) for 
specific functional forms of u, fi and fi. It is of interest to study the effects of the bath on 
these coherent states. In particular, in the absence of the bath, the coherent state evolves 
solitonically under the action of a Hamiltonian which is a linear combination of the gener- 
ators of the Lie group defining the coherent state ||. The coherent state is very special in 
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that it follows the exact classical equations of motion forever. However, since coherent states 
still display quantum interference it is incorrect to consider such states as truly classical. 
Much work has been done recently on the decoherence effects of an external bath and the 
suppression of quantum interference in the system |]]. More important is the study of the 
effects of the bath on quantum systems which, when isolated from the bath, do not possess 
coherent states. Does the interaction with the bath succeed in retarding the spread of an 
initial wave packet? ie. can one interpret the bath as effectively measuring the system and 
"collapsing" the wave function of the system? Does the interaction with the bath cause a 
deviation from the quantal equations of motion towards the classical equations of motion? 
As a first step towards answering these important questions one must solve for the quantum 
behavior of the system coupled to the bath. We do so in the quantum optical limit for a 
general system linearly coupled to the bath. 

The standard method used to solve equation (|1.1|) is to transform the master equation 
into a c-numbered partial differential equation for a particular quantum distribution func- 
tion (QDF). However, in certain important cases (degenerate parametric oscillator in the P 
representation), the resulting PDE can take the form of a Fokker-Planck equation with a 
negative diffusion matrix. Such cases do not correspond to the usual diffusive process. To 
use the machinery associated with diffusive processes one must make use of the generalized 
P-Representation [^]. This entails doubling the degrees of freedom and "choosing a gauge" 
to ensure that the diffusion matrix is positive and real. However, recent numerical work 
using this QDF has raised questions concerning it's validity ||. Of course, for system non- 
linearities higher than quadratic one has a c-numbered PDE of order three or more. Since 
this does not occur in this model we will not address this here. Resolving ( |1.2|) onto the x, p 



basis via a = (ux + ip) / \J27ilo we get 



AhuAp/-f = -u 2 {2n + l)[x, [x,p]] - (2n + l)[p, [p, p}} 

+iu)[p, px + xp] — iu[x, pp + pp] , (1.4) 



while resolving (|1.1| ) onto the Glauber-Cahill family of QDFs W(a,a*, s), which include the 



P (s=+l), Wigner (s=0) and Q (s=-l) quantum distribution functions we obtain 

W(a, a*, s) = - % - {[/i + 2a f 2 + a*u]d a * - [h + 2a* f 2 + au]d a - s[f 2 - f 2 ]d 2 aa * } W(a, a*, s) 
+ ^{d a a + d a *a* + {2n + l- s)d 2 aa ,}w{a,a*,s) . (1.5) 

From ( |1.4j ) we see that we now are effectively coupled to noise in both position and mo- 



mentum. We now convert equation ( |1 . 1|) into a Schroedinger equation on Super-Hilbert 
space. 

II. SUPER HILBERT SPACE AND LIE ALGEBRAIC SOLUTIONS 

We initially follow a technique developed by Barnett and Knight in association with 
the applications of Thermofield theory to quantum optics Q. In this theory, the density 
operator is transformed into a state vector in an expanded Hilbert space (Super-Hilbert 
space) and operations on operators become super-operators in the Super-Hilbert space. The 
transformation between Hilbert space and Super-Hilbert space is accomplished through the 
state 

I j> = X;ijv, ^ = 5:1^x1^ , (2.1) 

N N 

and through defining \p) = p\I). One can show a\I) = a'\I) and a^\I) = a\I) where the 
a's are the annihilation operators in Hilbert space 1 (the original space) and the a's are the 
annihilation operators in Hilbert space 2. Using (|2.1|) , equation (|1 . 1|) becomes a Schroedinger 
like equation on the Super-Hilbert space 

j t \p) = -iH\p) , (2-2) 



where 



- iH= -i [u [a^a - ~cJ~a] + f x [a - a f ] + A [a f - a] + f 2 [a 2 - S t2 ] + / 2 [a t2 - a 2 }} 

+ — (n + l)[2aa — a) a — a* a] + — n\2<v<v — a<v — a<v\ . (2-3) 
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The "Super-Hamiltonian" , H seems at first glance totally intractable, however, upon in- 
spection one can rewrite (|2.2|) as 

jU{t) = -iH(t) U(t) \p(t)} U(t)\p(0)} , (2.4) 

N 

hi{t)Ni , (2.5) 

8=0 

where the iVj are the generators of an N dimensional Lie algebra C We therefore have 
access to powerful Lie Algebraic techniques to solve (|2.4j ) for the propagator U(t). 

The algebra C can be identified as the two photon subgroup of Sp(6,R). This sub- 
group is a 15 dimensional, semisimple Lie algebra, possessing a five dimensional ideal X, 
{No, Ni, N2, N3, N4} and compact U{2) subgroup {N$, N$, Ni 3 , N u }. It also contains the 
single photon algebra and Weyl-Heisenberg algebra for both sets of oscillators. This 



algebra has been studied in detail by Gilmore and Yuan ||10|| . The description of the gener- 
ators Ni and coefficient functions hi(t) are given in Table | and Table [TT]. 

The standard Lie Algebraic method used to solve equations such as ( |2.4j ) is the Wei- 
Norman method [|ll|]- In this method one assumes a particular form for the propagator as 
an ordered product of exponentials of the generators 

N 

U{t) = Hexp(g i {t)N i ) . (2.6) 

i=0 



One then substitutes this ansatz into (|2.4j) to obtain differential equations relating g^ to 
hi. However, the resulting equations are highly non-linear and coupled. The complexity 
of these equations and the effort needed to obtain them grows rapidly with the dimension 
of the Lie group and proves prohibitive for semisimple Lie groups with dimensions greater 
than six. Even if one succeeds in solving the coupled ODEs for the g^s the global validity 
of the chosen ansatz is not guaranteed and must be tested through examination of the g^s 
|T2|| . Finally, the whole procedure is highly sensitive to the particular ordering of generators 
chosen in the ansatz for U(t). We will, instead, adopt a method discovered by Fernandez 
||13||, which can reduce the problem to the solution of a coupled set of linear ODEs and 



algebraic equations. Key to this approach is the existence of a suitably large proper ideal X. 
One effectively solves for the dynamics of the ideal and "lifts" this information via Baker- 
Campbell-Hausdorff (BCH) identities to determine the dynamics of the entire algebra. 

For any element of the proper ideal Ni G X (i G {0, .., 4}) we have U^IU G X where U 
is the propagator in fl2.4p. In particular 
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N t (t) = U-\t)NiU{t) = J2 u a(t) Nj , i G {0, ..,4} . (2.7) 

3=0 

It is easy to show that the time dependent coefficients Uy(t) obey 

N 

u = m-u , Hij{t) = £ M*)c& > u ij(t = °) = . ( 2 - 8 ) 

m=l 

where C^- are the structure functions of C. To solve for Uij(t) we must solve the linear set of 
coupled ordinary differential equations fl2.8p . By transforming to a new set of generators for 
L we can simplify the structure of (|2.8|) greatly. We will denote this new set by {Wj}. This 
change of basis decouples ( p.8|) into pairs of coupled linear first order differential equations 
and also greatly simplifies the BCH disentangling for the relations between the tty's and the 
coefficient functions gi(t) appearing in (|2.6|) . To obtain the time dependent coefficients gi(t) 



appearing in ( |2.6|) we use BCH disentangling identities to compute the functionals Fj(gi) 
appearing in 

U-\t)N t U(t) = £ Fj(gi)Nj . (2.9) 

3=0 



The F/s are algebraic functions of the ^'s. Equating the coefficients of Nj in (|2.7|) and (p.9|) 
for a particular Ni yields algebraic relations between the g^s and the Ui/s. To check the 
global validity of the particular ansatz chosen one must only determine when the algebraic 
relations become degenerate. 

Before we solve for the dynamics in the extended Super-Hilbert space we show first how 
to recover the dynamics of the density operator in the original Hilbert space. Resolving the 
density operator into number states gives 

P = Yl Pnm \n)(m\ , p nm = (n\p\m) . (2.10) 



We obtain the density state vector in the extended Hilbert space through 

\p) = p\I) = Pnm\n) x \m) = pnm\n,m) . (2.11) 

nm nm 

Thus, the original density operator may be recovered through 

p = (m, n\p) \n) (m\ . (2-12) 

nm 

Through similar steps one can also recover the components of the original density operator 
in the coherent state basis. We resolve unity in the coherent state basis through / = 
- / d 2 a \a)(a\ and using the identity 

oo 

Y: (P\n) \n) = \(3*) , (2.13) 

71=0 

we obtain 

\p) = -1 f d 2 a d 2 (3p(a* } (3) e -(l«l+l^l)/ 2 \ a , (3*) , (2.14) 

7T J 

where we have set |ot, (3) = \a) x \0) and p (a*,j3) = (a\p\{3). One can easily show 

(a\p\(3) = {(3\a\p) , (2.15) 
as expected. From Q2.15|) we can obtain the Q distribution function fl4[ 



Q(a*,a,t) = (a\p(t)\a) = (a*,a\p(t)) . (2.16) 

However, in the extended Hilbert space, \p(t)) = U(t)\p(0)). Using this with equations 
( P-14| ) and (|2.16|) we finally obtain 



Q(a*,a,t) = 

\ f d 2 ad 2 (3p (a*J)e-^ +W)/2 (a\a\U(t)\a,(3*} . (2.17) 



n 2 



In the examples treated in this paper the coherent state basis proves more useful than the 
number basis. Thus to calculate the Q function it will be necessary to compute 

(a*,a\U(t)\a,[3*) , (2.18) 



for various propagators U(t). In the systems treated, U will be at most quadratic in the 
annihilation and creation operators of the two Hilbert spaces. Writing 



\a*,P) = e-^ 2 -^ 2 e aat+/3 ^ t \0,0) 



(2.19) 



equation (|2.18f) becomes (0, 0|C/(t)|0, 0) where 



To evaluate the above we use BCH identities to normal order the operator U as 

fj(t\ = e R ii a i a j+ c i a i e Dij(a\a,j+5ij /2)+-qI ^Lijcuaj+lidi 



(2.20) 



(2.21) 



where the index i labels the two Hilbert spaces. The vacuum expectation value of U becomes 
trivial and yields 



(0,0\U(t)\0,0) = ex V ( Dll + D22 +r ] ) 



(2.22) 



III. SOLVING THE DYNAMICS 

In Table [III] we show the relation between Ni and Wi and in Table [TV] we give the 
commutation table for the ideal I. 

From the {Wi} commutation table we can construct the "Hamiltonian" Ti in 



Hi 



0, 0, 0, 0, 

-2h A , -h 8 /4-h 9 /2, -2h 6 , 0, 
-2h 3 , -2h 5 , -h 8 /4 + h 9 /2, 0, 



0. 

0. 



(3.1) 



2h 2 , 2h 7 , 0, hs/A-hg/2, 2h 6 , 

2h, 0, 2h 7 , 2h 5 , h s /4: + h 9 /2, 

We see that the degree of coupling within ( |2.8| ) is related to the number of nonzero entries 
in each row of Tiij. This is in turn related to the number of different W^s appearing in each 
row of the commutation table of the ideal X. From the initial condition Uij(t = 0) = 5^ and 
the structure of TCij we immediately obtain 



S 



u i0 = 5. 



iO 



U13 = U 23 = Uu = «24 = W 



Thus the set (12. 81) of ODE's for uu and u 2i can be re-cast as 



X x = zAXi + iB 
X 2 = iCX 2 + iD 



where 



Xo = 



( \ 

MlO Mil «12 



«20 «21 «22 



/ \ 

«33 «34 M31 «32 «30 



«43 M44 M41 M42 «40 



A = 



C 



B 



D 



^ 7^33 ^34 



7X21 7^44 



-z 7 /2 - w , 2/ 2 
y -2/2 , -i 7 /2 + w y 

-z 7 /2 -u 2f 2 ^ 
y -2/ 2 * 7 /2 + Wy 



w 10 » 



V 



n 2n 



^00 H 3 iu n H 3 iu 12 H30 ^ 



y 7^42^21 7^42^22 7^40 

From (|2.8|) , the initial value of Xi and X 2 are 



Xi(t = 0) 



x 2 (* = 0) 



'010^ 

yO 1 



^10 0^ 
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However, from Table |TT|, /14 = h$ = and thus B vanishes identically. Equation (|3.3|) then 
gives 

""10 = u 2o = Vt . (3.13) 

Note, however, that equation ( ^.3| ) must already be solved to correctly specify D in (|3.4|). 
The solution of equations (|3.3] , [3.4D together with (|3.2| ) completely determine the dynamical 
evolution of the 

To finally compute the propagator we must calculate the algebraic relations between the 
Qi(t) and Uij(t). We must evaluate the action of the group on each of the elements of the 
ideal, that is, we must calculate 

exp(- ft ad Wi) W j = e~ 9 ^ Wj e 9 ^ , (3. 14) 

for i G (0, ..,14) , j G (0, ..,4). This is achieved using Baker-Campbell-Hausdorff disen- 
tangling identities for Lie Groups. The relevant identities are given in Appendix 1. Tables 
|V] and [VT] summarize the action of the group on the elements of the ideal. The particular 
set of algebraic relations obtained will depend sensitively on the choice of ordering of the 
generators in the ansatz for U(t). We choose the ordering 

U(t) = f[.eM9i(t)Wi) , (3.15) 

where f]j is the product ordered from the left in the sequence 1, 2, 13, 5, 3, 11, 7, 8, 9, 14, 
10, 6, 12, 4, 0. With this ordering the six conditions ( |3.2| , |3.13|) give 

94 = 93 = 914 = , #io = -2#i2 , g u = ~2g 13 . (3.16) 

The very complicated disentangling of W u is thus avoided as g u = with this choice of 
ordering. Working through the BCH identities with a symbolic manipulator we finally arrive 
at the relations in Table |V1J| . 

The relations are analytic except when w n = or m 33 = 0. Thus the ansatz ( |3.15| ) is 
global if both uy 7^ and M33 7^ for all t. Except for go(t), the propagator is formally 
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determined. To solve for go(t) however, is not trivial. It cannot be determined through BCH 
identities as Wq commutes with every element of the group. Instead we must substitute 
( p,15| ) into the "Schroedinger" equation Ti, = iUU~ l and from this determine go. In p~3 



the resulting relation (just below equation (20) in [|13|]) is quite simple. However, for the 
much larger algebra Sp(6,r) the resulting relation can be quite complicated. To compute the 
above "Schroedinger" equation and obtain the general relation for g it is expedient to use 
the faithful matrix representation of Sp(6,r) given by Gilmore and Yuan [|K|. We will not 
give the general formula but will compute it case by case. 

IV. EXAMPLES 

In this section we apply the above methods to solve two model systems. The first is the 
simple harmonic oscillator. We set uj(t) = u and fx — f 2 — fi — f 2 — 0. We compute the 
matrix H,^ in (|2.8|) and construct the matrices A, B and C. We can easily solve the coupled 
sets of ODEs for X\ and X 2 to find the non-vanishing to be 

«ii = «i = u£* = nil = e (-^ 2 +^ , (4.1) 

u 31 = -2{2n + l)e~ iujt sinh — , u A2 = e 2iujt u 3l . (4.2) 

2 

The resulting non-zero g's are 

g 7 = (2n + l)e^ /2t sinh ^ , g 8 = 2 7 t , g 9 = -2iut . (4.3) 

The propagator on the extended Hilbert space is U(t) = e 97 w? e 98 Wa e" 9 e 90 °. This may 
be normal ordered in terms of SU(1,1) generators [[HJ, 

K + = a)a) , 
= aa , 
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Kq = a)a — a) a 



to give 



U(t) = exp(xK + ) exp(yK 3 + g 9 /2K ) exp(xKJ) exp(g W ) 

where 

= g 7 e 9s/4 - sinh g 8 /A 
X g 7 e9*/ 4 + cosh £ 8 /4 ' 

?/ = -2 In [# 7 e W4 + cosh # 8 /4] , (4.5) 

_ g 7 e 98//L + sinh g 8 /A 
Z ~ g 7 e3»/ 4 + cosh 8 /4 ' ^ ' ' 

To obtain g we compute 7i = iUU^ 1 to find go = 7/2. With U normal ordered and the 
identity 0, 

t^i LL a tv 



71 

Z=0 



it is a simple matter to compute ( 2.18|) , 



([3\a\U(t)\a,f3*} = 

1 



1 

—(5 /3) m a) exp 

K 



K 



{a-PxTiP-ax) 



(4.7) 



where k = 1 + n(l — e~ 7 *) and x = exp(— (7/2 + iu/)t). This result is well known and 
corresponds to the off-diagonal elements of the density matrix in the coherent state basis in 
the presence of a bath where p(t = 0) = |a)(/3|. 

For a second example we take u(t) = u> (l + Asm r/t), fi = ie luJdt D/2, fi = f* and 
f2 — f2 — 0. This corresponds to a variable frequency oscillator with a driving force. One 
again proceeds as before, however, the ODE's are slightly more complicated due to the 
forcing term. The non- vanishing g's are 

D/2(exp(iu p t) - exp(-7_t - e(cosr/t - 1))) 

9i = 77, ^ ; i : 7 > l 4 -°J 

7/2 — ibu + ^ smrjt 

= D/2(exp(-iuj p t) - exp(-7 + t + e(cosry - l)t)) 
7/2 + i5uj — - sin r/t 
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g 7 = (2n + l)e" 7 * /2 sinh ^ , (4.10) 
9s = 2lt , (4.11) 
g 9 = -2iu t + 2e(cos rjt - 1) , (4.12) 

where j± = 7/2 ± iu , 5u = u — u>d and e = iu A/rj. The propagator on the extended 
Hilbert-space is now U(t) = e 9lWl e 92 ™ 2 e 97 ™ 7 e 9sW * e " W9 e 9oWo . We again find g = 7/2. To 
compute the result analogous to ( |4.7|) we normal order the W7, Wg and Wg exponentials and 
insert a resolution of unity in the Super-Hilbert space between e 92W2 and e xK+ . To complete 
the calculation we perform the coherent state integral over the introduced partition of unity 

USing / £° e -Aa*a+Ba+Ca* = 1 The fin&1 regult j g 



0*,a\U(t)\a,fl*) 
— (a\/3) ((3\a) exp 

K 



—{gx + a - /3x)*(92 + (3 - ax) 

K 



(4.13) 



where 



In x = ~lt/2 + g 9 /2 = -jt/2 - i(J Q (t (cosr^t - 1)) 

V 



In the latter model U\\ = exp(7_t + e[cos rjt — 1]) and m 33 = exp(— j + t + e[cos rjt — 1]). 
Since \uu\ ^ and m 33 | ^ for all e and finite t the resulting propagator in both models is 
global. This method can thus treat complicated time dependencies. It can also be applied 
to very particular nonlinear systems |16| . 
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APPENDIX A: BCH IDENTITIES 



In this Appendix we derive the necessary BCH identities to disentangle the action of the 
propagator U (t) in ( |3.15| ) on the elements iVj of the ideal I. 
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From 17] we have the expansion 



e xA B e~ xA = B + x[A,B] + ^ [A, [A, B}} + 



(Al) 



where x is a c-number and A and B are elements of a Lie group. For the case where the 
commutator of A and B is a c-number we are left with the first two terms in the expansion 
(|Al~|). For the simplest non-trivial case where the commutation of A with B closes onto B 
itself, ie., [A, B] = mB one can show 



n x ad A 



B = Be Ti 



(A2) 



For the more complicated example of exp(— g u ad W\i) Wi the commutation does not close 



on the first iteration. Here we use Wilkox's method of parameter differentiation |18 
Letting 



xA r) „—xA 



G{x) = e xA Be 



(A3) 



with 



[A, B] = miB + m 2 C 



[A, C] = m 3 B + m A C , 



(A4) 



we can obtain 



An 

G\x) = ^ = [A, G(x)] 



(A5) 



We now assume that G(x) is of the form G(x) = a(x)A + b(x)B + c(x)C. With the initial 
condition G(x = 0) = B we see that a(x = 0) = c(x = 0) = and b(x = 0) = 1. 
Substituting this ansatz for G(x) into (|A5|) we obtain differential equations for a(x), b(x), 
and c(x). Immediately we get a(x) = and the coupled set 



b 



(A6) 



mi m 3 

For the particular example exp(— gi^ adWi^Wi the values for the rrii can be read off Table 
P"V| and yield ni\ = —1/2, irt2 = —1, m 3 = 2, = 1/2. Inserting these values into ( |Atj| ) and 
solving with the above initial conditions, gives 
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(A7) 



X 2 x ' X 
where \ = 2/ v7 ■ The results of the action of Wu on the elements of the ideal are given in 



Table IVTL 
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TABLES 

TABLE I. We give the realisation of the two-photon subalgebra of Sp(6,r) in terms of creation 
and annihilation operators a) ,a) , a, a. 



N 


iVi 


N 2 N 3 


N A N 5 N 6 N 7 


N s Ng 


Nio Nn N12 Ni 3 


JV 14 


1 


a 


a a* 


a' aa a a 


a' a aa' 


a* 2 a 1 " 2 a* a* a* a + 1/2 


2*3 + 1/2 






TABLE II. 


Table of coefficients hi, 


appearing in 


the 'Super-Hamiltonian' (|2,5|). 




ho 


hi 


h2 hz 


hi hs ha hj 


hs hg hio 


hii hi2 his 


/214 


7 /2 


h 


-f* n 


-fl i-y(n + 1) / 2 -/ 2 * 


/* 


—/a «7n a; — 2772 — 17/2 — a) - 


- ijn — i'y/2 
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TABLE III. We give the relation between the new generators W% and the old generators N{. 



Also given 


are the new coefficient funct 


ions appearing in (pi 


) and the corresponding two-photon 


operator. 








Wi 


N; 


hi 


Operator 


Wo 


A^o 


7 /2 


1 


Wi 


iVi -N 4 


fx 


a — a* 


w 2 


N 2 - N 3 


-n 


a — a) 


w 3 


N 1 +N 4 





a + a* 


w 4 


N 2 + N 3 





a + a* 


w 5 


N 6 - N n 


h 


a 2 -at 2 


W 6 


N 7 - N w 


-fl 


at 2 -at 2 


w 7 


N 5 + N 12 - {N 13 + N M ) 


i 7 (2n + l)/2 


aa + at&t — (at a + atg + 1) 


w 8 


(N 5 - N X2 )jA 


2i 7 


(aa - atat)/4 


w 9 


(N 13 - N 1A )/2 


2w 


(at a -at a) /2 


Wio 


N 8 





a' a 


W u 


N 9 





a) a 


w 12 


N 7 + N 10 





a 2 + at 2 




N 6 + iVn 





a 2 + at 2 


W14 


2N 5 + N 12 - (N 13 + N 14 )/2 





2aa + a f a f - (at a + a) a + l)/2 



TABLE IV. The commutation table for the ideal J = {Wo, Wi, W 2 , W 3 , W 4 } with the whole 
two-photon group £. The entries are [Wj,Wj]. 



W t \Wj 


Wi 


w 2 


W 3 


Wa 


W 5 


W 6 


W 7 


W 8 


Wg 






W12 


Wis 




Wi 











2W 





2W 2 





Wi/4 


Wi/2 


W 4 





2W 4 





Wi/2 + W 3 


w 2 








2W 





2Wi 








W 2 /4 


-W2/2 





n 3 





2W 3 


W 2 /2 + W 4 


w 3 





-2W 











-2W 4 


-2Wi 


-W 3 /4 


Ws/2 


H 2 





-2W 2 





-2Wi ~ W3/2 


W4. 


-2W 











-2W 3 





-2W 2 


-W4/4 


-W4/2 





Hi 





-2Wi 


-2W 2 - Wa/2 
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TABLE V. Table of the action of the group element Wj on an element Wj of the ideal J, i.e. 
each entry shows exp(— giWi)Wj exp(gjWj). 



Wi\Wj 


Wi 


W 2 


W 3 


W 4 


Wi 


VFi 


w 2 


W 3 


W 4 - 20i 


w 2 


Wl 


w 2 


W 3 - 2g 2 


W 4 


w 3 


Wi 


W 2 + 2 53 


W 3 


W 4 


w 4 


Wi + 254 


w 2 


w 3 


W 4 


w 5 


Wi 


W 2 + 2g 5 W 1 


w 3 


W 4 " 2 55 ^3 


w 6 


Wi + 2(7 6 W2 


w 2 


W 3 - 2g 6 W 4 


w 4 


w 7 


Wi 


w 2 


W 3 - 2g 7 W 1 


W 4 - 2g 7 W 2 


W 8 




W 2 e 98 / 4 


W 3 e -98 / 4 




VVg 




W„ o — 99/2 


W„ o99/2 

KK 3 6 a ' 


W. o — 99/2 


Wio 


Wi + S10W4 


W 2 


W 3 + 510^2 


W 4 


Wii 


Wi 


W2+5HW3 


w 3 


W 4 + g xx W x 


W12 


Wi + 2012 W 4 


W 2 


W 3 - 2g 12 W 2 


W 4 


Wi 3 


Wi 


w 2 + 2 5l3 w 3 


w 3 


W 4 - 2<7i3 Wi 



TABLE VI. Table of the action of W14 on the elements of the ideal I. We have set 



tan 


Hi = 7 x l 2 and [i 2 = gu/a where a = 2/7 1 / 2 . 


See the Appendix for more details. 




W u 


Wi 




0(^/2 cos(^i - fi 2 )Wi + sin ^2^3) 


w 2 




a(\/2cos(^i - n 2 )W 2 + sin/^Wi) 


w 3 




a(\/2cos(^i + /U 2 )W3 - 2 sin ^2^1) 


w 4 




a(\/2cos(^i + /U 2 )W4 - 2sin fi 2 W 2 ) 



19 



TABLE VII. Algebraic relations between the coefficients gt and the functions uij. 

gi = -U40/2 g 2 = -«3o/2 53 = 

94 = g 5 = U2i/(2un) 96 = ui2/(2un) 

97 = -«3i/(2un) 5s = 21n(nn/n 33 ) 59 = ln(unu 33 ) 

gw = -2g 12 gu = -2g 13 g 12 = - [u 32 - u 31 u 12 /uu] /(4u 33 ) 

513 = - [«41 + «21«3l/«ll] / (4nn) 
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